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the exchange transitions which occur at small 
internuclear distances, the rate coefficients for the 
three surfaces agree quite well. For the direct 
transitions, there are significant differences for the 
pure rotational transitions from j=0 to 2 and from j=l 
to j = 3 in which there is no change in vibration. For 
higher j the differences tend to disappear, suggesting 
that the rotational angular momentum can couple to the 
orbital angular momentum to overcome the centrifugal 
barrier so that the excitation is driven by a torque 
operating at small separations of H and H 2 where the 
three potential energy surfaces are similar. 

There are of course substantial questions as to the 
validity of semi-classical calculations. We have also 
performed complete numerically exact quantum mechanical 
calculations for the processes in which vJ changes from 
00 to 01, 02 and 03 and from 01 to 00, 02 and 03. The 
calculations are computationally intensive and have been 
carried out only for the DMBE surface and for 
temperatures up to 1000K. Table 2 is a comparison of 
semi-classical and quantum calculations for the rate 
coefficients of the (0,0) — » (0,2) and (0.1) to (0,3) 
transitions in crn^s -1 . It appears that the semi- 
classical calculations are useful, indeed accurate, down 
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to temperatures of about twice the threshold temperature 
of the reaction which is 509K for (0,0 )— > (0,2) and 944K 
for (0,1) — > (0,3) . 

A first draft of a paper on the quantal 
calculations has been written and is attached. 

We have been studying the problem of the 
photodissociation of NH . Our preliminary study has 
identified the mechanisms and has definitely established 
that the lifetime is of the order of 10 5 s but it also 
demonstrates that a more detailed analysis is needed. 
Dr. M. A'Hearn has provided us with his data on the 
fluorescent population of the NH rotational and fine- 
structure levels from which we should be able to predict 
accurate photodissociation lifetimes. 

We have also under study an investigation of the 
distribution rate of C 2 but have not so far identified 
the mechanism responsible. For CO, we are carrying out 
explicit calculations. The literature value is based on 
inadequate laboratory data. 

A review of H 3 + in terrestrial and extraterrestiral 
environments was prepared for a volume of Advances in 
Atomic, Molecular and Optical Physics in honor of Sir 
David Bates. A copy is attached. 
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Table la. State to state rate coefficients for H4H 2 at 1000K with the LSTH potential. 
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Table 2. State to state rate coefficients for H+H 2 
with the DMBE potential 


temperature 


(0,0) 

to (0,2) 

(0,1) 

to (0,3) 


This 

paper 

Sun and Dalgarno 

This paper 

Sun and Dalgarno 

200 

2.0 ( 

-13) 

1.9 (-12) 

1.7 (-16) 

l.H-13) 

400 

7.1 ( 

-12) 

1 . 5 ( -11 ) 

9.6 (-14) 

2.1(-12) 

600 

2 .3 ( 

-11) 

3 .4 (-11) 

9.5 (-13) 

6.7 (-12) 

800 

5.4 { 

-11) 

5.4 (-11) 

4.0 (-12) 

1.3 (-11) 

1000 

7.6 ( 

-11) 

7.3 (-11) 

7.8 (-12) 

2 .0 (-11) 
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Abstract 

The rate coefficients for the excitation from the j — 0 to the j = 2 
rotational level of H 2 and from the j = 1 to the j — 3 rotational level 
by impacts with hydrogen atoms are calculated for temperatures be- 
tween 30K and 1000K using a fully- converged complete close-coupling 
method. Rate coefficients are also obtained for the reactive process 
of ortho-para conversion where j is excited from j = 0 to j — 1,3 
and from j = 1 to j = 2. The interference between the direct and 
reactive channels is taken into account as is the geometric phase re- 
sulting from the adiabatic separation ot electronic and nuclear motion 
that generates the potential energy surface. Convenient analytical 
representations of the rate coefficients are presented. 
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1 Introduction 


The distribution of molecular hydrogen between its ortho and para forms 
is an important aspect of the physics and chemistry of the atmosphere of 
the Jovian planets (Massie and Hunten 1982, Conrath and Gierasch 1984), 
and of interstellar molecular clouds (Dalgarno, Black and Weisheit 1973, 
Spitzer, Cochran and Hirshfeld 1974, Jura 1975, Black and Dalgarno 1977, 
Takavanagi, Sakimoto and Onda 1987, Pineau des Forets, Flower and Mc- 
Carroll 1991, Abgrall et al. 1992). Depending on the physical environments, 
conversions between ortho and para forms may be effected by reactive colli- 
sions between H and IF in which the rotational quantum number j changes 
by an odd number. Rotational excitations in which j changes by an even 
number in non-reactive and reactive collisions of H with H 2 are important to 
the distribution of the relative populations of the rotational levels of H 2 and 
to the cooling efficiency of hydrogen rich gas. 

Because of its fundamental significance to the theory of chemical reac- 
tions, the reactive scattering of H by H 2 has received considerable attention 
and powerful methods have been developed for and applied to the calculation 
of the cross sections (Schatz and Kuppermann 1976, Mladenovic et al. 1988, 
Zhang and Miller 1988, 1989, Park and Light 1989, Walker and Light 1989, 
Manolopoulos and Wyatt 1989, 1990, Launay and Dourneuf 1989, Miller and 
Zhang 1991). Most of the effort has been devoted to an exploration of the 
dynamical resonances that may or may not have been detected experimen- 
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tally at energies near leV (iSieh and Valentini 1988, 1990, Rliner, Adelman 
and Zare 1991), but some calculations of reactive scattering at lower temper- 
atures have been carried out (Park and Light 1989, Takayanagi and Masaki 
1991) and Mandy and Martin (1991, 1992, 1993) have used a semi-classical 
procedure to calculate rate coefficients for all the possible reactive and non* 
reactive rovibrational transitions in H 2 for temperatures down to 600K. 

No calculations of rotational transitions at lower temperatures have been 
reported since the calculations of Green and Truhlar (19/9) who did close- 
coupling calculations oi non-reactive rotational excitation processes within 
the rigid rotor approximation, neglecting any interference between reac- 
tive and non-reactive channels. They obtained rate coefficients considerably 
smaller than the results of earlier close-coupling calculations of Allison and 
Dalgarno (1967), Wolken, Miller and Karplus (1972), Chu and Dalgarno 
(1975), and McGuire and Kruger (1975). 

In many astrophysical applications, rate coefficients determined by Elitzur 
and Watson (1978) by interpolation and extrapolation of the early calcula- 
tions (Allison and Dalgarno 1967, Nishimura 1968) have been used. 

It is now practical to carry out numerically exact calculations of the 
rotational excitation cross sections at low velocities which do not use the 
rigid rotor approximation and which include all the important reactive and 
non-reactive channels. We employ here the generalized Newton variational 
principle GNVP (Schwenke et al 1988, Sun et al. 1989) to calculate the cross 
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sections for the j = 0 — * 1,2,3 and j = 1 — > 0,2,3 transitions of H 2 in its 
lowest vibrational state with H at total energies below 0.9eV. We include all 
the contributing states, both open and closed, and obtain the non-reactive 
and reactive S matrix elements in the same calculation. The non-reactive 
calculations of Green and Truhlar (1979) were based on the LSTH potential 
energy surface (Truhlar and Horowitz 1978, 1979). We make used of the 
modified DMBE version ( Varandas et al. 1987) which is designed to be more 
accurate at large separations of H and H 2 . 

2 Theory 

In the Arthurs-Dalgarno (1960) representation, the cross section for for ro- 
tational transition of a diatom from rotational level j to rotational level j l 
by atom impact at a total energy E is given by 

<7 = (1) 

where kj is the wave number of the relative motion of the colliding particles 
and P(E,j — » j') is given in terms of the S-matrix element S J for the J th 
partial wave by the expression 

J=0 
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The S-matrix for collisions of H and H 2 is given by 


\S J (E,jl - j'V ) | 2 = ~2\S J n {E,il - jT)\ 2 , 

3\Sj{E,jl - jT)\\ 

1 \Sl(E,jl — j'l')\ 2 t 


\Si(EJl - j'l') + Sj{E,jl — ;T)| 2 , j,j' even; 
|5^(£,j7 - m - 5/(£,jZ - j'7')| 2 + 


j,j' odd; 
j even, j' odd; 
j odd, j ; even; 

( 3 ) 


where 5^ and 5/ are the non-reactive and reactive S-matrix elements re- 
spectively. Eq. (3) was obtained from formalism of Miller (1969) except that 
here we have reversed the sign oi S? in front of S/ to compensate for the 
geometric phase arising from the adiabatic separation of the electronic and 
nuclear motion (Mead and Truhlar 1979, Mead 1980, Lepetit and Ivupper- 
mann 1990). The permutation symmetry for three identical identical nuclei 
leads to an irreducible representation, which reduces the size of the problem 
to one third (Miller 1969, Schwenke et al 1989). 

The rate coefficient is defined as 

i(r,j - ;') = - a (4) 

where ks is Boltzmann’s constant, T is the temperature, fi is the reduced 
mass of the colliding system, H-H2, Ej is the rotational energy of level j of 
H2, and Eo is the larger value of Ej and Ej » . 


3 Calculations 

We calculated the S-matrix up to a total energy of 0.9eV, varying the basis 
function size to ensure numerical convergence to better than 5%. As the 
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energy increases, more partial- waves contribute, requiring more extensive 
calculations. 

Reactive processes only occur at low values of the angular momentum 
J at which the projectile may overcome the centrifugal barrier and enter 
the short-range reaction region. At high J. the reactive channels may be 
neglected and the calculation simplifies to a non-reactive problem which re- 
quires much less computational effort. Near E = 0.8eV, for example, the 
reactive transition amplitudes become very small beyond J = 15 but the 
non-reactive ones do not until J = 40. According to Eq. (3), the integral 
cross section has more contributions from high J S-matrix elements due to 
the weighting factor 2 J + 1 and is less sensitive to the reactive processes that 
occur at low J. This is not true for differential cross sections where the low 
J contributions may dominate at particular scattering angles. In practice, 
we found also that below 1Q00K it is a valid approximation to neglect the 
closed ro- vibrational channels and the reactive channels for the para-para 
and ortho-ortho transitions. 

Fig. I shows the calculated cross sections for the j = 0 — ► 1,2,3 and 
j = 1 — 0,2,3 transitions, from which we obtained rate coefficients with an 
error of less than 15%. We used cubic-spline fitting for the calculated P(E ), 
or log P(E) in the lower energy region, where P(E) grows exponentially. 

Table I lists the rate coefficients k(j — ► j*) for the non-reactive rotational 
excitations j = 0 — ■ 2 and j = 1 — > 3. The rate coefficients of the inverse 
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transitions can be obtained according to 

•^ '1 E' -E- f 3 -7 even and j' odd, 

k{j' — j) = ~t~— ~ reip( )fc(j — j") x ( 1/3 j odd and j' even, 

^ 1 ^ l 1 otherwise. 

( 5 ) 

The energies of the rotational levels of H 2 are (Ej — Eo)/kB == 1 i OAT, 509 AT , 
and 1014A' for j = 1,2,3 respectively. 

The values for the 0 — 2 and 1 — '3 transitions are much larger than the 
rate coefficients calculated by Green and Truhlar (1979) which are also listed 
in Table I. Although the calculations of Green and Truhlar (1979) used the 
rigid rotor approximation and neglected closed channels and reactive chan- 
nels, there seems little question that the large discrepancies arise primarily 
from the different potential energy surfaces employed in the two calculations. 
Given that they used a semi-classical method, the calculations of Mandy and 
Martin (1993), who also adopted the LSTH surface but did not use the rigid 
rotor approximation or neglect the reactive channels, are in acceptably close 
agreement with the results of Green and Truhlar (1979). Their values of 
k(0 — 2) and k(l — 3) at 1000K are included in Table 1. 

Further evidence is provided by unpublished preliminary semi-classical 
calculations by S. Lepp (private communication 1993) who obtains much 
larger values of k(Q — * 2) and k( 1 — * 3) for the DMBE surface than for the 
LSTH surface. 

The surfaces differ mostly in the long range region, to which the non- 
reactive cross sections are sensitive at low temperatures. It is claimed by 
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Varandas et al. (1987) that their potential is superior to the LSTH potential 
at large separations of H and H 2 . If so, the rate coefficients we have presented 
in Table 1 are preferable to those of Green and Truhlar (1979) and Mandy 
and Martin (1993). By chance, they are of the same order of magnitude as 
the early calculations of Allison and Dalgarno (1967). 

The rate coefficients may be conveniently represented as functions of tem- 
perature by the expressions 


k ( 0 ~ 2) = exp(— 509/T) x 10” 13 cxn 3 s” 1 

k( 1 — 3) = exp( — 844/T) x 10“ 13 cm 3 s“ 1 


(0.86T -f 36), 30A" > T < 100 A, 

1.3T, 100A" > T < 1000A, 

(0.29T + 5.9), 30 K >T< 100 A, 

(0.47T - 14), 100A >T< 1000 A, 


to within the accuracy of the calculated rate coefficients. 

The rate coefficients for the reactive excitation and de-excitation pro- 
cesses are presented in Table 2. In contrast to the non-reactive excitation 
process, the values agree well at least at the common temperature of 1000K 
with the semi-classical results of Mandy and Martin (1993). At 1000K we 
obtain k ( 0 — > 1) = 7.8 x 10 _13 cm 3 s _1 , k ( 0 — + 3) = 3.9 x 10~ 13 cm 3 s -1 
and k( 1 — + 2) = 2.1 x 10~ 13 cm 3 s"" 1 , whereas they calculated respectively 
6.4 x lO" 13 ^ 3 ^ 1 , 4.7 x 10' 13 cm 3 s _1 and 1.9 x 10 _13 cm 3 s” 1 . The different 
potential energy surfaces used are quite similar at shorter distances of ap- 
proach where reaction occurs. The agreement supports the utility of the semi- 
classical method above 600K (Mandy and Martin 1993). Takayanagi and 
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Masaki (1991) have used a simplified ^-conserving coupled states method. 
With the LSTH surface they obtain k(l — > • 0) = 1.9 x 10~ 18 cm 3 s 1 and 
4.4 x 10 -21 cm 3 s -1 at 200IC and 100K, while we obtain 9.8 x 10 _19 cm 3 s* 1 
and 5.4 x 10 _22 cm 3 s _1 . At low temperatures, the J, conserving method is 
probably inadequate though some of the differences are due to the potentials. 
Experimental data are available for thermal rate coefficients between 300K 
and 444 K. Schultz and LeRov (1965) give an analytic fit to the data, repro- 
duced in the review of Schofield (1967). We have complete results only for 
initially populated j = 0 and j = 1 levels so that no comparison is possible. 

The uncertainties in the calculated rate coefficients may be as large as 
15% at the lowest and largest temperatures. For the transitions from j = 0 
to j =■ 1,3 and j = 1 to j = 2, the analytical expressions: 

k(0 -r l) = exp[(5.30 — 460/T) 2 — 21.2]10 -13 cm 3 s -1 , (6) 

k(0 — 3) = expi(6.36 — 373/T) 2 — 34.5ll0 -13 cm 3 s -1 , (7) 

k{ 1 — 2) = exp[(5.35 - 454/T) 2 - 23.1]10- 13 cmV 1 . (8) 

provide fits to our results to within the accuracy of the calculated rate coef- 
ficients. 
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Figure captions 

Fig.l: the cross sections for the transition from j = 0 to j = 2 and from 
j = 1 to j = 3. 


Fig. 2: the cross sections for the transition from j = 0 to j = 1,3 and from 
j = 1 to j = 2. 
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Table 1 . Rate coefficients for the 0 — 2 and 1 — * 3 rotational excitation 


t r an sition of H; by H impact* 


T(K) | 

0 — 2 i 

1 — 3 

i i 

a b c 

a b c 

i 

1 30 

2.7( - 19) 

9.2(-25) 

i 40 

2.1(-17) 

1.2(-21) 

50 

3 . 0 ( - 16) 

9.4(-20) : 

60 

1 . S ( - 1 5 ) 

1 -8(- 18 ) 

70 : 

6.6( - 15) 

1 -5(- 1 7 ) 

: 80 

1 . 8 ( - 1 4 ) 

7 -5(- 17) 

90 

4.0(-14) 

2.7(-16) 

100 

7.7( - 14 ) l-0(- 15) 

! 7.7(-16) 7.5(- 17 ) 

150 ! 

6.2(-13) 

' 1.9(-14) 

200 

i 1 . 9 ( - 1 2 ) 

1.K-13) 

300 

7 . 2 ( - 1 2 ) 9.7(- 14 ) 

7.4(-13) 4.4(-14) 

400 

l-5(-ll) 

2. 1(- 12 ) 

500 

2.4(-ll) 3.6(- 13) 

4.1(-12) 2.0(-13) 

600 

3.4(-ll) 

6.7(-12) 

700 

4.4(*11) 

9.7(- 12) 

800 

5.4(-ll) 

1 -3(- 11) 

900 

i 6.4(- 11) 

1.6(-11) 

1000 

7.3(-ll) 2.1(-12) 9.4(- 12 ) 

2.0(- 11) l-2(-12) 1.2(-12) 

i units of cm 3 sec x ; number in parent] 

ieses is the power of 10; 


(a) this work, (b) Green and Truhlar (1979), (c) Mandy and Martin (1993). 


15 


Table 2. Rate coefficients for the 0 — *1,3 and 1 — *2 rotational excitation 


transition of H 2 by H impact* 


T(K) 

0 - 1 

0 — 3 

1 — 2 

100 

1.0 (-22) 

1.0 (-24) 

1.8 (-23) 

! 150 

8.5 (-21) 

5.1 (-22) 

1.9 (-21) 

200 

4.2 (-19) 

5.7 (-20) 

1.1 (-19) 

300 

9.3 (-17) 

2.2 (-17) 

2.3 (-17) 

400 

2.1 (-15) 

6.1 (-16) 

5.2 (-16) 

| 500 

1.5 (-14) 

5.0 (-15) 

3.7 (-15) 

600 

5.6 (-14) 

2.1 (-14) 

1.4 (-14) 

700 

1.4 (-13) 

6.0 (-14) 

3.7 (-14) | 

o 

o 

oo 

2.9 (-13) 

1.3 (-13) 

7.6 (-14) j 

900 

5.1 (-13) 

2.4 (-13) 

1.3 (-13) j 

1000 

7.8 (-13) 

3.9 (-13) 

2.1 (-13) I 


in units of cm 3 sec number in parentheses is the power of 10. 
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Emission lines of H 3 + in the 2u 2 <2) -> 0 overtone band of H 3 + are present in the energetic electrons. They are transformed into H 3 + ions by reactions (I) (4). The H 3 + ions arc 
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understanding of the absorption spectrum of H 3 + , carried out an experiment using an infrared magnetospheric ions. The mechanism that accelerates the ions presents a challenging problem in 
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absorption technique, in which the recombination of H 3 + ions in the v =0 level was measured. He plasma physics. 




[ 3 + + e -> H 2 + H (5) obtained a rate coefficient of 1 . 8 x 10 " 7 cm^ s * 1 at 300K. The interpretation of his data was subject 
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The rate coefficients *12 and *13 are known functions of the temperature (Smith, Adams and 
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